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ABSTRACT 

We present the results of a search for mid-infrared signs of star formation activity in the 1.1 mm 
sources in the Bolocam Galactic Plane Survey (BGPS). We have correlated the BGPS catalog with 
available mid-IR Galactic plane catalogs based on the Spitzer Space Telescope GLIMPSE legacy survey 
and the Midcourse Space Experiment {MSX) Galactic plane survey. We find that 44% (3,712 of 8,358) 
of the BGPS sources contain at least one mid-IR source, including 2,457 of 5,067 (49%) within the 
area where all surveys overlap (10° < I < 65°). Accounting for chance alignments between the BGPS 
and mid-IR sources, we conservatively estimate that 20% of the BPGS sources within the area where 
all surveys overlap show signs of active star formation. We separate the BGPS sources into four 
groups based on their probability of star formation activity. Extended Green Objects (EGOs) and 
Red MSX Sources (RMS) make up the highest probability group, while the lowest probability group 
is comprised of "starless" BGPS sources which were not matched to any mid-IR sources. The mean 
1.1 mm flux of each group increases with increasing probability of active star formation. We also find 
that the "starless" BGPS sources are the most compact, while the sources with the highest probability 
of star formation activity are on average more extended with large skirts of emission. A subsample of 
280 BGPS sources with known distances demonstrates that mass and mean H2 column density also 
increase with probability of star formation activity. 

Subject headings: star: formation — infrared: stars — dust — ISM: clouds 



1. INTRODUCTION 

The formation of massive stars remains one of the 
prominent puzzles in astronomy (see review by McKee & 
Ostriker 2007). A current goal in high-mass star forma- 
tion studies is to determine an evolutionary sequence and 
lifetimes similar to the Class system in low-mass star for- 
mation (e.g. Shu et al. 1987). While the evolutionary se- 
quence is still open to considerable debate, one proposed 
sequence is summarized in recent reviews by Churchwell 
(2002) and Zinnecker & Yorke (2007). This evolutionary 
sequence begins with prestellar cores or clumps: gravi- 
tationally bound over-densities within a molecular cloud 
which show signs of inward motion but have not yet be- 
gun to form a protostar. Since prestellar cores are lacking 
any internal heating, their temperatures are typically 10- 
20 K and their spectral energy distributions peak in the 
far-infrared. 
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Infrared dark clouds (IRDCs), seen in absorption 
against a diffuse mid-infrared background, are consid- 
ered the cold precursors to star clusters (Rathborne ct 
al. 2006), and thousands have been cataloged from the 
MSX and GLIMPSE surveys (Egan et al. 1998; Carey et 
al. 1998; Simon et al. 2006; Peretto & Fuller 2009). While 
IRDCs may be the sites of the earliest stages of massive 
star formation, the clumps and cores within IRDCs can 
host massive young stellar objects (MYSOs) in various 
stages of evolution (e.g. Chambers et al. 2009; Rathborne 
et al. 2005, 2006, 2007, 2008, 2010). While a prestellar 
core is only seen as an IRDC if it is located at the near 
distance and there is enough background IR emission to 
absorb, IRDCs are seen in emission at far-infrared, sub- 
millimeter and millimeter wavelengths, allowing for de- 
tection and identification of prestellar cores throughout 
the Galaxy. 

The second stage in the proposed evolutionary se- 
quence of Churchwell (2002) is hot cores: compact (di- 
ameter < 0.1 pc), dense {n > 10^ cm"'^), and warm 
(T> lOOK) molecular cloud cores with high molecular 
line brightness temperatures (Kurtz et al. 2000) which 
are a result of internal heating from a massive protostar 
(Churchwell 2002). These hot cores will not contain any 
radio free-free emission because the infall rate will be 
high enough for the infalling material to absorb the UV 
flux from the central protostar. 

Once the accretion begins to taper off, a small but 
detectable HII region will form. Thus, the third evo- 
lutionary stage is thought to include Hyper-Compact 
HII (HCHII) regions which will rapidly expand to Ultra- 
Compact HII (UCHII) regions. This stage can be char- 
acterized by the presence of radio emission (Churchwell 
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2002). The final stages are compact and classical HII 
regions. The HII region has now expanded and begun 
disrupting the parent molecular cloud while revealing the 
embedded stellar population in the optical and infrared. 

This evolutionary sequence based on the size of HII 
regions may be too simplistic. Recent simulations of the 
gravitational collapse of a massive molecular cloud sug- 
gest that the size and morphology of HII regions can 
vary drastically during the formation of a massive star, 
with the HII region expanding and contracting repeat- 
edly through the HCHII and HII classifications (Peters, 
et al. 2010a; 2010b). These simulations suggest that the 
size of an HII region will only depend on age late in the 
lifetime of the HII region once accretion has halted. 

Other evolutionary sequences have been proposed 
based on the presence of various types of masers, ra- 
dio emission, thermal emission from warm dust, as well 
as any combination or complete lack of these star forma- 
tion tracers (e.g. De Buizer et al. 2005; Minier et al. 2005; 
EUingsen et al. 2007; Longmore et al. 2007; Purcell et 
al. 2009; Breen et al. 2010). These studies have included, 
at most, on the order of a hundred sources, and while 
detailed studies of small samples are crucial for under- 
standing the results of large-scale surveys. Galaxy-wide 
samples are needed in order to solidify an evolutionary 
sequence. 

Recent and ongoing (sub) millimeter wavelength con- 
tinuum surveys of the Galactic plane are providing a 
novel inventory of star formation sites within the Galaxy 
at unprecendented sensitivity and resolution (BGPS, 
Aguirre et al. 2011; ATLASGAL, Schuller et al. 2009; 
JPS; Di Francesco et al. 2008). These surveys are de- 
tecting thousands of sites of massive star formation and 
providing a statistically significant sample with which to 
begin studying star formation on a Galactic scale. How- 
ever, these surveys alone cannot provide direct evidence 
of a source's evolutionary state and whether it is actively 
forming stars. 

Fortunately, a wealth of other multiwavelength surveys 
that can provide evidence of a source's evolutionary state 
are already available or currently being observed. For 
example, the Spitzer Space Telescope GLIMPSE survey 
(Benjamin et al. 2003; Churchwell et al. 2009) at 3.6 to 
8.0 /um can provide information regarding the popula- 
tion of embedded young stellar objects (YSOs) found 
throughout the Galactic plane via mid-infrared colors 
(e.g. Robitaille et al. 2008) and the presence of outflows 
from massive young stellar objects traced by extended 
4.5 /im emission (e.g. Cyganowski et al. 2008). Similarly, 
the Red MSX Source survey (RMS; Hoare et al. 2004) 
has identified a reliable sample of massive YSOs by fol- 
lowing up color-selected MSX YSO candidates. The 
Spitzer MIPS Galactic Plane Survey (MIPSGAL; Carey 
et al. 2009) has mapped the emission from warm dust in 
the Galactic plane at 24 and 70 ^m and can also identify 
embedded protostars. The on-going Herschel Infrared 
Galactic Plane Survey (Hi-Gal; Molinari et al. 2010) will 
map the same region as the Spitzer surveys at 60-600 
/im and will identify YSOs at wavelengths spanning the 
peaks of their spectral energy distributions. Together 
with the shorter wavelengths from Spitzer^ Hi-Gal will 
provide accurate protostellar luminosities (e.g. Elia et 
al. 2010). At longer wavelengths, the Co-Ordinated Ra- 
dio 'N' Infrared Survey for High-mass Star Formation 



(CORNISH; Purcell et al. 2008) will provide a complete 
catalog of UCHII regions. To gain the most insight from 
the copious Galactic plane surveys, we must begin cor- 
relating their data sets. 

In this paper, we present the results of a comparison of 
the Bolocam Galactic Plane Survey (BGPS) and avail- 
able. Galaxy-wide mid-IR based YSO catalogs with the 
goal of characterizing the star formation activity within 
the BGPS sources. This work represents the first step in 
determining the evolutionary state of the BGPS sources. 
Further comparison to other Galactic plane surveys, such 
as MIPSGAL and CORNISH, is required and wiU be 
the subject of future papers. Section [5] presents the de- 
tails of the various data sets included. The basic results 
are presented in Section [3l including the cross- matching 
method, a discussion of the prevalence of chance align- 
ments, and the basic cross-matching statistics. Section 
[4] discusses trends in various millimeter properties as a 
function of star formation activity, and Section [S] pro- 
vides a summary of the work. 

2. DATA SETS 

2.1. The Bolocam 1.1 mm Galactic Plane Survey 

The Bolocam Galactic Plane Survey^" (BGPS; Aguirre 
et al. 2011) has surveyed 170 square degrees of the North- 
ern Galactic Plane in 1.1 mm continuum emission us- 
ing the facility long-wavelength camera, Bolocam (Glenn 
et al. 2003; Haig et al. 2004), at the Cahech Submil- 
limeter Observatory (CSO) on Mauna Kea. The sur- 
vey consists of two distinct portions: a blind survey of 
the inner Galaxy spanning —10.5° < I < 70.0° where 
|6| < 0.5° everywhere except for 1.0° sections in I cen- 
tered at ^ = 3°, 15°, 30°, and 31° where |6| < 1.5°, and 
a targeted outer Galaxy portion including Cygnus-X 
(70° < £ < 90.5°, |5| < 1.5°), the Perseus Arm (i - 
111°), IC1396 {I -99°, 5-3.5°), the W3/4/5 region 
[I - 135°,5 - 0.5°), and the Gemini OBI Molecular 
Cloud {I - 190°,6 - 0.5°). 

Bolocam is a hexagonal array of 144 bolometers with a 
filter centered at 268 GHz with an effective band center 
of 271.1 GHz (Aguirre et al. 2011) and a bandwidth of 
45 GHz. The filter excludes the C0(2-l) emission line. 
Each individual bolometer has a Gaussian beam with a 
FWHM of 31", and the instrument has an instantaneous 
field of view of 7.5' due to the separation of the bolome- 
ters in the focal plane of the instrument. For more infor- 
mation regarding Bolocam, see Glenn et al. (2003), and 
Haig et al. (2004). 

A reduction pipeline was created specifically for the 
BGPS survey, and is described in Aguirre et al. (2011). 
In order to retain the extended structure in the final 
maps, an iterative sky cleaning method was employed. 
At each iteration a model of the sky emission is up- 
dated and subtracted from the overall data timestream 
to reveal a time stream containing only the astrophys- 
ical signal. The astrophysical signal is then mapped 
with 7.2" X 7.2" pixels and the effective beam size for the 
BGPS is 33". 

Observations were obtained with scans along both I 
and h using the raster scan mode and a scan speed of 
120"s~^. Sensitivity to extended structure up to 5.9' an- 
gular size was retained by not utilizing chopping. The full 

See |http: / / milkyway.colorado.edu/bgps/ 1 



3 



flux density of objects larger than 5.9' is not recovered, 
but long filaments with at least one dimension smaller 
than 5.9' are still detected (Aguirre et al. 2011). The flux 
density was calibrated based on observations of planets. 
Aguirre et al. (2011) compared BGPS with 1.2 mm data 
acquired with MAMBO on the IRAM 30 m (Rathborne 
et al. 2006; Motte et al. 2003, 2007) and SIMBA on the 
SEST 15 m (Matthews et al. 2009) and found that the 
vl.O BGPS data were systematically lower and require 
a mean scaling factor of 1.5±0.15 to align the BGPS 
fluxes with the previous MAMBO and SIMBA data sets. 
Here we present the vl.O data with the scaling factor of 
1.5±0.15 applied. 

Millimeter continuum sources were extracted from the 
iterativly mapped images using the BGPS source extrac- 
tion software, Bolocat (Rosolowsky et al. 2010). Bolo- 
cat employs a seeded watershed method (Soille 1999) 
which first identifies all pixels above a signal-to-noise ra- 
tio of 2.0, identifies contiguous groups of pixels above 
this signal-to-noise ratio, culls any groups which contain 
fewer than a beam's worth of pixels (22 7.2" pixels), and 
then expands each remaining high signal-to-noise region 
through a nearest-neighbor algorithm to include all adja- 
cent pixels above a signal-to- noise ratio of 1.0. Extending 
the region to include pixels down to 1 u is based on the 
reasoning that marginal emission next to regions of sig- 
nificant emission is likely real. The software then breaks 
each region up into smaller substructures based on the 
contrast between local maxima. Through the nearest 
neighbor algorithm each source is assigned a group of 
contiguous pixels and a source number. Bolocat then 
produces a label map, in which every pixel assigned to 
a source is assigned the source number. Source proper- 
ties are then computed via emission-weighted moments. 
All flux dependent properties presented in this work are 
based on the integrated flux density which is simply the 
sum of all pixels assigned to a BGPS source. This flux 
density avoids contamination from nearby sources that 
can occur with aperture photometry. For more details 
regarding source extraction and parameter estimation, 
see Rosolowsky et al. (2010). 

In the generation of the Bolocam catalog (Rosolowsky 
et al., 2010) we evaluated both the completeness limit 
(false negatives as a function of signal-to-noise ratio) and 
the spurious detection incidence (false positives as a func- 
tion of catalog parameters) . The completeness study was 
presented as part of the original catalog manuscript. The 
catalog parameters were specifically chosen to minimize 
false detections as well as false negatives. To evaluate 
the false-positive incidence, we inserted an array of 2600 
false sources into the data time stream for a blank field 
of observation. The time stream was for the 1=111 field 
which had been cleaned of signal as a byproduct of the 
iterative mapping procedure (Aguirre et al. 2011). The 
injected sources had a range of brightnesses and the re- 
covery fraction as a function of brightness established 
the completeness limit. The algorithm recovered false 
detections - all sources in the final catalog were the in- 
jected sources. While this may be surprising, the cata- 
log's relatively high effective completeness limit of 5a is 
a byproduct of eliminating false detections. Taking this 
as ± 1 over the 4 deg^ area of the field leads to an up- 
per limit of ^ 40 false detections in the survey, or 0.5%. 
We note that the injected sources were point sources and 



extended source recovery tests found that the BGPS cat- 
alog software often broke up extended emission into mul- 
tiple objects, but it did not lead to the false detection of 
sources not associated with emission. Based on this test, 
we believe the BGPS catalog contains only a very small 
number of false sources. 

The BGPS catalog is >99% complete at the 5a level. 
While a varies across the survey due to the varying num- 
ber of observations and observing conditions, the sur- 
vey as a whole is >98% complete at the 0.4 Jy level 
(see Figure 9 in Rosolowsky et al. 2010). Assuming a 
dust temperature of 20 K (Dunham et al. 2010), this 
completeness level translates into a minimum mass of 
5 Mq, 21 Mq, 130 Mq, 524 Mq, and 750 M© at 1 
kpc, 2 kpc, 5 kpc, 10 kpc, and 12 kpc, respectively. 
However, since the BGPS filters out diffuse emission on 
scales larger than 5.9' (Aguirre et al. 2011), the types 
of sources the BGPS is sensitive to depends on distance 
(see Section 7.2 of Dunham et al. 2010 and Section 6 of 
Rosolowsky et al. 2010). At the nearest distances (up 
to approximately 2 kpc), the BGPS will be sensitive to 
cores (mass M= 0.5 — 5 M0, volume-averaged density 
n= 10^ — 10^ cm~'^, and physical radius R= 0.01 — 0.1 pc; 
Bergin & Tafalla 2007), while at intermediate distances 
the BGPS wiU be sensitive to clumps (M= 50 - 500 Mq, 
n= 10^ - 10"* cm-^, R= 0.15 - 1.5 pc; Bergin & Tafalla 
2007), and at the furthest distances BGPS will detect 
entire molecular clouds (M= 10^ - 10* Mq, n= 50 - 500 
cni-^ R= 1 - 7.5 pc; Bergin & Tafalla 2007). 

The BGPS images and catalog have been released to 
the public and are hosted by the Infrared Processing and 
Analysis Center via the NASA/IPAC Infrared Science 
Archive" (IPAC). 

2.2. GLIMPSE Red Source Catalog 

In addition, we made use of the GLIMPSE survey, a 
Spitzer /IRAC Legacy survey of the Galactic mid-plane 
(Benjamin et al. 2003; ChurchweU et al. 2009) at 3.6, 
4.5, 5.8, and 8.0 ^im. The GLIMPSE I survey covers 
10° < \£\ < 65° and |6| < 1° and the GLIMPSE II 
survey fills in the region for \£\ < 10° with \b\ < 1° 
for 1^1 > 5°, |6| < 1.5° for 2°< \i\ < 5°, and |6| < 2° for 
\e\ < 2°. The total area of the GLIMPSE I and II surveys 
is 274 deg^, and the overlap region with the BGPS survey 
is -10.5°< i < 65°, |6| < 0.5° and |6| < 1.0° at i =3°, 
15°, 30°, and 31°. The data products include both highly 
reliable Catalogs, and less reliable but more complete 
Archives. 

Robitaille et al. (2008; hereafter R08) compiled a 
highly reliable census of 18,949 GLIMPSE point sources 
with excess emission at mid-infrared wavelengths, of 
which 11,649 were classified as candidate young stellar 
objects (YSOs), and 7,300 as candidate asymptotic giant 
branch (AGB) stars. The AGB classification was further 
separated into two categories consisting of "standard" 
AGB stars and "extreme" AGB stars. The "extreme" 
AGB stars have high mass-loss rates and therefore sig- 
nificant circumstellar dust. In this paper, we refer to 
the "standard" and "extreme" AGB stars as sAGB and 
xAGB stars, respectively, and both types together as sim- 
ply AGB stars. The separation of the sources into YSOs 

" See |http://irsa.ipac.caltech.eclu/data/BOLOCAM_GPS/| 
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and AGB stars was based on IRAC and MIPS color and 
magnitude selection criteria, and was only approximate. 
In reality, 50 — 70% of the sources are likely to be YSOs, 
and 30 — 50 % are likely to be AGB stars. 

The R08 census was designed to maximize reliability. 
As the sensitivity to point sources decreases in areas of 
bright diffuse emission - which are common at 8 /im - 
the census used stringent magnitude selection criteria so 
as not to be affected by this. As a result, the ROB census 
may be missing fainter YSOs that are still detected by 
GLIMPSE. Therefore, for the current work, we compile 
a list of red GLIMPSE sources that were missed in the 
R08 census, with the caveat that this complementary list 
is therefore less reliable, and is background rather than 
flux limited, especially at 8 /im. 

We constructed this complementary sample as follows. 
We first increased the reliability of the GLIMPSE catalog 
fluxes in a similar way to ROB. First, each position in the 
GLIMPSE survey was observed at least twice, and we 
therefore removed fluxes from the Catalog that depend 
only on a single detection. We required the fractional 
rms of the individual detections and the fractional error 
in the combined flux to be less than 15%, otherwise the 
flux was removed. We then rejected sources that were 
not detected at least at 4.5 and 8.0 ^m, and finally we 
selected the red sources by requiring [4.5] — [8.0] > 0.75, 
a less stringent color selection than that used for ROB 
([4.5] - [8.0] > 1.0). This results in 64,859 sources being 
selected, excluding the ROB sources, which we refer to as 
additional GLIMPSE sources. 

Robitaille et al. (200B) estimated that the red source 
catalog is 65.7% complete. Additionally, they explored 
the range of luminosities of YSOs detectable as a func- 
tion of distance. The sensitivity of the GLIMPSE survey 
results in a lower flux limit, and Robitaille et al. imposed 
an upper flux limit in order to exclude point sources with 
unreliable fluxes due to saturation. They used the mod- 
els of Robitaille et al. (2006) to calculate average [4.5] 
and [B.O] magnitudes for Stage I and II YSOs as a func- 
tion of luminosity, and calculated the range of distances 
at which a YSO would be included in their Red Source 
Catalog. A 1 Lq YSO within 1 kpc will be included in the 
red source catalog, while a 10 Lq YSO will be included 
if it is between 0.2 and 2 kpc. A 10^ Lq YSO between 
2 and 10 kpc will be included in the R08 catalog, and a 
10^ L0 YSO between 4 and 11 kpc will also be included. 
The most luminous YSOs will only be included at the 
farthest distances due to the upper flux limit imposed 
to avoid saturation, while the lowest luminosity sources 
will only be found at the nearest distances due to the 
sensitivity limits of the GLIMPSE survey (see Figure 19 
and Section 4.5 of Robitaille et al. 2008). 

2.3. Extended Green Objects Catalog 

The GLIMPSE Catalogs and Archives arc point source 
catalogs, and therefore exclude any extended YSOs. In 
particular, massive embedded YSOs driving outflows can 
in some cases be seen as extended sources at 4.5 /im, 
as this band contains H2 emission lines as well as CO 
band heads (Reach et al. 2006, Smith et al. 2006). These 
lines can be excited by shocks, which are expected to 
be present where outflows interact with the interstel- 
lar medium. Cyganowski et al. (200B) compiled a list 



of 304 such objects, dubbed Extended Green Objects"'^^ 
(EGOs), from the GLIMPSE I survey The EGOs are 
likely tracing outflows from massive YSOs since they 
show a strong correlation with IRDCs and Class II 
methanol masers (Cyganowski et al. 2009), both of which 
are associated with early stages of massive star formation 
(Cyganowski et al. 200B). The overlap area with the 
BGPS survey is 10°< t < 65°, 1&| < 0.5°and \b\ < 1.0° 
at i =15°, 30°, and 31°. 

The EGO catalog is certainly not a complete census 
of YSOs with shocked outflows. The requirement of ex- 
tended 4.5 /im emission likely resulted in the preferential 
inclusion of YSOs at relatively small distances, as the 
angular extent decreases with distance. Also, as noted 
in Cyganowski et al. (2008), EGOs near bright sources or 
in areas with signiflcant PAH emission and diffuse green 
emission extended on scales less than 10" are likely ab- 
sent from the catalog. 

Thus, the EGO catalog may include more massive, 
higher luminosity, nearby sources which were excluded 
from the ROB catalog because of the saturation limits. 

2.4. Red MSX Source Catalog 

The Red MSX Source (RMS) survey^^ (Hoare et 
al. 2004; Urquhart et al. 2008b) was designed to follow-up 
and conflrm massive young stellar objects (MYSO) can- 
didates selected from the Midcourse Space Experiment 
(MSX) Galactic plane survey (Price et al. 2001). The 
Galactic plane survey included all longitudes from to 
360°, and latitudes within \b\ < 5° and consisted of imag- 
ing at B.28, 12.13, 14.65, and 21.3 /im with a resolution 
of IB" at aU wavelengths (Price et al. 2001). A total of 
1,992 MYSO candidates were identifled using color selec- 
tion criteria. Low- and intermediate-mass YSOs were ex- 
cluded from the catalog. All MYSO candidates were sub- 
sequently followed-up to conflrm their nature (Mottram 
et al. 2007; Urquhart et al. 2007a, 2007b, 2008a, 2009a, 
2009b) and were classified into the following categories: 
"Evolved Star", "PN", "OH/IR Star", "Young/Old? 
Star", " HII Regions", "HII/YSO", and "YSO". In this 
work we include the "PN", "OH/IR Star", and "Evolved 
Star" classifications in a general "Evolved Star" category. 
The "YSO/HH" classification is assigned when both an 
HII Region and YSO are present or if a definitive classi- 
fication is not possible, while the "Young/old" classifica- 
tion is assigned to sources which look like evolved stars 
but have strong CO emission along their line of sight re- 
sulting in an uncertain classification. The BGPS overlaps 
the RMS survey where £ > 10°. 

Unlike the ROB sample, the RMS selection criteria did 
not include an upper-limit on fluxes, and as such the 
highest luminosity sources at all distances will be in- 
cluded in the catalog. YSOs with L'^ 10'^ Lq will be 
included to 3 kpc, while 10* Lq YSOs will be included 
to a distance of 10 kpc, and 10^ Lq YSOs will be in- 
cluded across the entire Galaxy (see Figure 2 of Urquhart 
et al. 2008b) . The luminosity range of the RMS survey 
overlaps with the range in the ROB red source catalog, 
and fills in the region of luminosity space excluded from 

The 'green' adjective comes from the fact that these objects 
appear as green in IRAC three color images, since these often use 
4.5 fim for the green channel 

See http://www.ast.leeds.ac.uk/RMS] 
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Fig. 1. — Left: The gray scale image shows the label map from a small portion of the BGPS. Right; The corresponding BGPS image 
shown in inverted grayscale ranging from 0.2 to -0.1 Jy/beam. In both images, the black ellipses mark the positions and sizes of the 
extracted BGPS sources as determined from emission weighted moments. Crosses and diamonds represent GLIMPSE R08 red sources 
which are and are not associated with a BGPS source, respectively. Plus signs and circles mark additional GLIMPSE sources which were 
matched and not matched with a BGPS source, respectively. Note the BGPS source near {i,b) ~ (6.23°, —0.39°). The ellipse encloses 
many pixels which do not include millimeter emission, thus demonstrating the benefit of using the label map for matching rather than the 
mean source properties which are shown here as the black ellipses. 



the R08 catalog because of the saturation limitations. 

3. RESULTS 

3.1. Cross-Matching Method 

The data sets were cross-matched based on spatial co- 
incidence. All infrared sources located within a BGPS 
source as defined by the label maps (see are consid- 
ered to be associated with the BGPS source. The label 
maps produced as a part of the BGPS source extraction 
(Rosolowsky et al. 2010) represent the full extent of each 
1.1 mm continuum emission source down to the la level 
of the original image. We require that the catalog posi- 
tion of each mid-IR source fall within a pixel assigned to 
a BGPS source. The EGO sources subtend more than 
a single BGPS pixel; however we find that ignoring the 
extended emission and simply using the catalog position 
corresponding to the peak of the 4.5 /im emission does 
not adversely affect the matching statistics for the EGO 
catalog. Figure [1] shows the label map and correspond- 
ing image centered near I 6° with the positions of the 
matched (crosses, plus signs) and unmatched (diamonds, 
circles) R08 and additional GLIMPSE sources, respec- 
tively. The catalog parameters for each BGPS source 
(major and minor axes and position angle) are displayed 
as the black ellipses. The BGPS sources clearly do not 
fill the black ellipses, highlighting the benefit of using the 
label maps rather than the BGPS mean source properties 
when cross-matching catalogs. 

3.2. Chance Alignments 

In matching various catalogs of mid-IR sources to the 
BGPS sources, we need to take into account the possi- 
bility of chance alignment of mid-IR sources with BGPS 
clumps. Since we do not have distance estimates to the 
mid-IR sources or most BGPS sources, we cannot deter- 
mine which mid-IR sources that lie along the same line of 



sight as a BGPS source are in fact chance alignments. In- 
stead, we can estimate the fraction of chance alignments 
statistically. 

In order to do this, we first consider the fraction of 
the BGPS area that contains mm emission and assume 
that the mid-IR sources are evenly distributed across the 
survey area. In this case, the fraction of area covered 
with mm emission corresponds to the fraction of mid-IR 
sources that are chance alignments with the mm con- 
tinuum emission. We calculate this fraction for the re- 
gion that is common to all surveys: 10° < i < 65°, and 
|6| < 0.5° with extensions to |6| < 1.0° as described in 
i i2.ll We find that 2.1% of the BGPS area contains source 
emission, and therefore 2.1% of the mid-IR sources will 
overlap mm continuum emission by chance. 

In reality, the distribution of the mid-IR sources is not 
even across the survey areas, but is concentrated towards 
the Galactic mid-plane, and also decreases with increas- 
ing longitude (Figure [2]) ■ Furthermore, even a popula- 
tion of mid-IR sources unrelated to star formation could 
display a similar non-uniform spatial distribution. To 
obtain a more realistic estimate of chance alignments, 
we attempt to preserve the large-scale distribution of the 
mid-IR sources by duplicating the catalogs and random- 
izing each source's coordinates within a circle of radius 
1° (in cases where mid-IR sources fell out of their original 
survey area, their positions were re-randomized) . In this 
way, any population unrelated to star formation would 
preserve its large-scale distribution and would not show a 
difference in the number of chance alignments, while any 
star forming population that was originally preferentially 
aligned with BGPS emission on arcsecond to arcminute 
scales would lose such a correlation. Figure [3] shows 
the distributions of Galactic longitude (left) and latitude 
(right) for the GLIMPSE red sources with £> -10° for 
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Fig. 2. — The distribution in Galactic longitude (left) and Galactic latitude (right) of the BGPS sources (solid hue), GLIMPSE Red 
Sources with I > —10° (dashed), and all RMS sources (dash-dot). The EGO catalog is not plotted due to the small number of sources 
(84). The sources are clearly concentrated toward the Galactic mid-plane (6 ~ 0°), and also show structure in Galactic longitude. 



the original (solid lines) and randomized (dashed lines) 
catalogs. Overall, the trends in the distributions of the 
original catalog are preserved in the randomized catalog 
although the peak in sources toward 6 ~ 0° is slightly 
decreased. To investigate the effect of the randomization 
radius used, we repeated the experiment with a radius of 
0.1° and found that the fraction of chance alignments in- 
creases slightly. All GLIMPSE R08 source types (xAGB, 
sAGB, and YSO) taken together have 8.1% chance align- 
ments with a 1° randomization. If we use only a 0.1° 
randomization, we expect 10.7% chance alignments. The 
choice of radius for the randomization does not severely 
alter the chance alignment results. 

Multiple randomized catalogs were created for each 
mid-IR data set in the interest of obtaining an error bar 
on the number of chance alignments. For the GLIMPSE 
catalogs, we randomized the catalogs 10 times, while for 
the EGO and RMS catalogs we randomized the catalogs 
100 times to avoid issues due to small-number statistics. 
Column 5 of Table 1 presents the percent of matched 
mid-IR source types that are expected to be chance align- 
ments, and column 9 shows the percentage of BGPS 
sources expected to be matched due to chance align- 
ments with each mid-IR source type. Overall, around 
7% of mid-IR sources would be expected to match the 
BGPS emission by chance alignments. 

3.3. Matching Statistics of the mid-IR sources 

Columns 2-5 of Table [T] present a summary of the 
matching statistics for each mid-IR catalog, including 
both the number and fraction of mid-IR sources that are 
matched to BGPS emission. The EGOs are most often 
matched to BGPS emission, with 90% of EGO sources 
matched (more than ten times that expected from ran- 
dom chance alignment). In comparison, the RMS sources 
have a lower overall matched fraction (68%). However, 
if we consider the different types of RMS sources sep- 
arately, it becomes apparent that the types of RMS 
sources associated with star formation (e.g. HII Re- 
gions, HII/YSO, and YSO) have a much higher fraction 
of matches with BGPS sources. Table [T] lists the match- 
ing statistics for the RMS sources by category. The 
Young/old star, HII Region, HII/YSO, and YSO cate- 
gories have matched fractions approximately 8 — 10 times 
the fraction expected by chance alignments (80-90%). In 



contrast, the fraction of RMS evolved stars matched with 
BGPS emission is perfectly consistent with that expected 
from chance alignments. Of the 133 RMS sources clas- 
sified as AGB stars within the "Evolved stars" category, 
only 12 were matched with BGPS sources. 

The GLIMPSE R08 YSOs tend to be less often as- 
sociated with mm emission (33% - or 4 times that ex- 
pected from chance alignments) . This low fraction could 
be due to several reasons. First, Stage II YSOs, which are 
present in the GLIMPSE ROB catalog, no longer have an 
envelope and are unlikely to possess the density and an- 
gular size required for a detection in the BGPS (Dunham 
et al. 2010, §7.2). Secondly, the separation of YSOs and 
AGB stars in the GLIMPSE Red Source catalog was only 
approximate, and refining the magnitude and color selec- 
tion criteria could result in a higher fraction of matched 
YSOs. The fraction of AGB stars from the GLIMPSE 
R08 catalog matched to mm emission is slightly above 
the level expected for chance alignments. This could be 
due to the approximate YSO/ AGB separation which re- 
sults in contamination from genuine YSOs in the sources 
classified as AGB stars, or it could be indicative of some 
AGB sources being present in the BGPS catalog. In- 
deed, the dust around some AGB and post-AGB stars 
is detected at millimeter wavelengths (e.g. Buemi et 
al. 2007; Dehaes et al. 2007; Ladjal et al. 2010) and it 
is possible that some BGPS sources are actually evolved 
stars. The RMS team has performed follow up observa- 
tions and has identified the sources in their catalog which 
are AGB stars. Although the number of RMS "Evolved 
stars" matched with a BGPS source is consistent with 
the number expected due to chance alignments, we can 
use the matching statistics to estimate the number of 
GLIMPSE R08 AGB stars expected to be matched with 
a BGPS source. If we assume the matching statistics for 
the RMS AGB stars within the "Evolved stars" category 
(12 of 133 were matched with a BGPS source, 9.0%) also 
apply to the R08 AGB stars, then we expect 296 R08 
AGB stars to be matched with BGPS sources. We find 
a total of 494 R08 AGB stars are matched with BGPS 
sources, 198 more than expected based on the matching 
statistics of the RMS AGB stars. If the RMS statistics 
are representative of the ROB AGB stars as well, then 
roughly 60% of the matched ROB AGB stars are likely to 
be AGB stars and 40% are likely to be YSOs. Further 
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Fig. 3. — Distributions of Galactic longitude (left) and latitude (right) of GLIMPSE red sources with I > —10° for the original (solid 
lines) and randomized (dashed lines) catalogs. Randomizing by 1° has very little effect on the distribution of I, but does slightly wash out 
the peak in the distribution of h. Overall, the trend of more sources located near the Galactic mid-plane (6 ~ 0°) is maintained. 



observations are required to determine the nature of the 
AGB sources which are matched with BGPS sources. A 
future paper will further explore the matching statistics 
in order to assess the validity of the magnitude and color 
cuts made to define the YSO GLIMPSE Red Source list 
presented in R08. 

Finally, the additional GLIMPSE source catalog con- 
tains 21% of sources that are matched with mm emis- 
sion (almost three times the level expected from chance 
alignments), indicating that it does contain a fraction of 
genuinely young stars. 

3.4. Star Formation Activity in BGPS Sources 

In columns 6 — 9 of Table [l] we present the detailed 
numbers/fractions of BGPS sources containing each type 
of mid-IR source. We assume the following mid-IR 
source types indicate star formation within the matched 
BGPS source: GLIMPSE R08 YSOs, sAGB, and xAGB 
stars; GLIMPSE additional sources; EGOs; young RMS 
source types including "Young/old star", "HII Region", 
"HII/YSO", and "YSO". As previously mentioned (sec- 
tion [331), '^'2 use the R08 AGB catalogs as an indicator 
of star formation because of the approximate separation 
of the YSOs and AGB sources and the resulting contam- 
ination of the ROB AGB catalogs by YSOs (see section 
12. 2p . Further observations would be required to deter- 
mine if R08 AGB sources matched with BGPS sources 
are YSOs or AGB stars, and we note that some AGB 
stars matched with BGPS sources will actually be AGB 
stars. 

A total of 3,712 BGPS sources were identified as con- 
taining signs of star formation via the alignment of a 
mid-IR source, with the caveat that some of the matched 
R08 AGB stars are actually AGB stars and not indica- 
tive of star formation. This equates to 44% of the entire 
BGPS catalog, although this is a lower-limit to the true 
fraction since the mid-IR catalogs do not overlap com- 
pletely with the BGPS catalog. Only the RMS catalog 
overlaps the BGPS at ^ > 65°, and unlike the GLIMPSE 
survey it will only identify the brightest mid-IR sources. 
If we consider only the region where all surveys overlap 
(10° < I < 65°), we find a total of 2,457 BPGS sources 
containing a mid-IR source out of 5,067 BGPS sources 
in the inner Galaxy (48%). 

In order to calculate the true number and fraction of 



BGPS sources that show signs of active star formation, 
we must account for the occurance of chance alignments. 
As previously discussed, the number of BGPS sources 
likely to have chance alignments were calculated for each 
mid-IR source catalog (the product of columns 6 and 9 
in Table [IJ. We simply subtract the number of BGPS 
sources identified as chance alignments from the number 
of BGPS sources which were matched with each indi- 
vidual mid-IR source type, and find that 1,472 of the 
8,358 (17.6%) BGPS sources in the entire catalog are 
reliable matches. If we restrict this calculation to the 
5,067 BGPS sources in the area where all surveys over- 
lap (10°< e <65°), we estimate that 1,035 (20.4%) BGPS 
sources are reliably matched with a mid-IR source. These 
percentages are conservative estimates of the percentage 
of BGPS sources that are actively forming stars and are 
driven so low due to the large number of BGPS sources 
expected to be identified via chance alignments with the 
additional GLIMPSE sources (1,573 for the entire cata- 
log; 976 for the overlap area only). 

Since we are unable to identify which mid-IR sources 
are chance alignments with our current data set, we con- 
sider all matches in the subsequent discussion. Table [5] 
presents the number of mid-IR sources within each BGPS 
source. Column 1 gives the BGPS source name. Columns 
2 and 3 give the R.A. and DEC of the geometric centroid, 
while Columns 4 and 5 give the Galactic coordinates. 
The angular radius of each BGPS source is given in Col- 
umn 6. The number of matched sources from the EGO, 
RMS, GLIMPSE R08, and additional GLIMPSE cata- 
logs are given in Columns 7-10, respectively. The full 
table is included in the online journal only. 

We define a Confidence value for each type which is 
the percent of matched BGPS sources which are not due 
to chance alignments. This is defined as the number of 
BGPS sources with a matching mid-IR source minus the 
number expected from chance alignments, divided by the 
number of BGPS sources with a matching mid-IR source, 
and is given in column 10 of Table [1] 

This confidence value allows us to estimate the prob- 
ability that a match with a particular mid-IR source in- 
dicates ongoing star formation in the BGPS source. The 
calculated confidence values range from 26%, for BGPS 
sources matched with RMS evolved stars, to roughly 
90%, for the BGPS sources matched to EGO and young 



8 




Fig. 4. — Number of sources as a function of BGPS source flux for the different probability groups of BGPS sources. Only the 5067 
BGPS sources in the region where all catalogs overlap are included. Note the vertical scale is different for each probability group. The 
range of confidence values for each probability group are given in the upper-right corner of each panel. 



RMS sources. The confidence is around 70% for the 
GLIMPSE R08 YSOs, around 50% for the additional 
GLIMPSE sources, and around 50% and 40% for the 
GLIMPSE R08 xAGB and sAGB stars respectively. 

4. DISCUSSION 

In this section, we consider only the 5,067 BGPS 
sources located in the region where all catalogs overlap 
(10°< I < 65°). In order to compare the properties of 
the BGPS sources as a function of the star- formation 
activity they contain, we define a classification scheme 
based on the confidence parameter listed in Table [TJ We 
separate the BGPS sources into four groups representing 
the probability of star formation activity: 

3. BGPS sources containing at least one EGO 
or young RMS (Young/old star, HII Region, 
HII/YSO, or YSO) source (corresponding to a con- 
fidence of l80%) 

2. BGPS sources containing no EGO or young RMS 
source, but at least one GLIMPSE R08 YSO (cor- 
responding to a confidence of 70%) 

1. BGPS sources containing no EGO, young RMS 
source, or GLIMPSE R08 YSO, but at least one 
GLIMPSE ROB AGB star or additional GLIMPSE 
source (corresponding to a confidence of 40-50%). 
The inclusion of the R08 sources classified as AGB 
stars may seem counter-intuitive, but they are con- 
taminated by genuine YSOs and are therefore to 
some extent a star formation indicator. 

0. The remaining BGPS sources, which show no sign 
of star formation in the above catalogs. We refer 
to these sources as "starless" . These may in fact be 
forming stars that are too faint to be detected, or 



too young and too deeply embedded to be visible 
in the mid-IR. 

In this classification scheme, a BGPS source is assigned 
to a group based on the highest confidence mid-IR source 
type it was matched with. For example, a BGPS source 
matched with an EGO, R08 YSO, and 5 additional 
GLIMPSE sources would be assigned to group 3 based 
on the presence of the EGO source, while a BGPS source 
without the EGO, but with a R08 YSO and 5 additional 
GLIMPSE sources would be assigned to group 2 based 
on the presence of the R08 YSO. 

Probability Groups through 3 are comprised of 2,610, 
1,324, 793, and 340 BGPS sources, respectively. The 
trend of fewer BGPS sources with increasing confidence 
parameter refiects the smaller number of mid-IR sources 
in the catalogs required for each higher confidence cate- 
gory. We note that the actual probability of active star 
formation for each BGPS source is unknown. The prob- 
ability we refer to throughout this paper is a measure of 
the probability that the mid-IR source is a real associa- 
tion rather than a chance alignment. 

4.1. Observed Properties 

Figure m shows the number of BGPS sources as a func- 
tion of the fiux density integrated over the entire source 
for the BGPS sources in each probability group. BGPS 
sources in Group 3 with the highest probability of star 
formation activity have the highest mean fiux, 8.0 Jy, 
with a standard deviation about the mean of 18.0 Jy. The 
mean fiux decreases with star formation activity proba- 
bility, with a mean of 0.71 Jy and standard deviation of 
1.4 Jy for the "starless" sources. The fiux distributions 
are asymmetric with tails extending to large fluxes. For 
example, the distribution of "starless" sources has a high 
flux tail extending to approximately 30 Jy. 
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Fig. 5. — Number of sources as a function of the fraction of 1.1 
mm flux density which is due to free-free emission based on 6 cm 
radio continuum observations toward RMS sources (Urquhart et 
al. 2009c) as discussed in section 14.11 The mean fraction of flux 
density attributed to free-free emission is 0.039, and the median is 
0.013. Overall, free-free emission is not the cause of the trend of 
increasing flux density seen with increasing probabiUty of star for- 
mation, although for some BGPS sources the free-free contribution 
can be significant (> 50%). 



The 1.1 mm BGPS flux densities could be contami- 
nated by free-free emission from ionized gas if the BGPS 
source is coincident with an HII region or Planetary 
Nebula, for example. The trend of increasing flux den- 
sity with probability of star formation could be due to 
an increase in the contamination from free-free with in- 
creasing probability of star formation. To quantify the 
free-free contamination we have compared the BGPS 
flux densities with integrated intensities from 6 cm radio 
continuum observations toward RMS sources (Urquhart 
et al. 2009c). We have identified 223 radio continuum 
sources from Urquhart et al. which coincide with 159 
unique BGPS sources. In cases where multiple radio con- 
tinuum sources coincide with a BGPS source, we have 
simply taken the sum of the individual 6 cm fluxes for 
comparison to the 1.1 mm flux. Assuming the 6 cm emis- 
sion is on the optically thin side of the free-free spectrum, 
we scale the 6 cm emission to 1.1 mm via Si,{S, 1.1mm) — 
S'i/(6cm)(l.lmm/60mm)'''^ since S,j{bee — free) oc v~'^'^ . 
If the 6 cm emission is on the optically thick side of the 
free- free spectrum, then 5*1/(5,1.1 mm) underestimates 
the free-free flux density, as the spectrum increases as v'^ 
on the optically thick side (Mezger & Henderson 1967; 
Tielens 2005). 

Figure [5] shows the distribution of the ratio of the free- 
free flux density at 1.1 mm to the 1.1 mm BGPS flux 
density. The dotted line marks the mean fractional con- 
tribution from free-free of 0.039 ± 0.066, and the dashed 
line marks the median of 0.013. The minimum and max- 
imum free-free contributions are 0.0002 and 0.54, respec- 
tively, and both values are for BGPS sources in the "star- 
less" group. The RMS 6 cm observations were targeted 
toward potential HII regions with the goal of further clas- 
sifying the RMS sources; some RMS sources were classi- 
fied as Planetary Nebulae which we did not consider as 



a positive sign of star formation and could therefore fall 
into our "starless" probability group. The vast major- 
ity of the BGPS sources including 6 cm emission fall in 
probability group 3, with a mean free-free contribution 
of 0.033±0.046, median of 0.012 and a maximum of 0.23. 
For all but four of the 159 BGPS sources with 6 cm ra- 
dio observations, contamination from free-free emission 
does not exceed 20%, while for the remaining four BGPS 
sources the free-free is contributing up to approximately 
half of the 1.1 mm flux density. 

Without radio observations for all BGPS sources, we 
are unable to correct for the free-free contamination. To 
assess if the mean level of free-free contamination we see 
towards the 159 BGPS sources could be the cause of the 
increasing flux density with probability of star forma- 
tion, we calculate what fraction of the group 1-3 flux 
would need to be attributed to free-free emission for the 
mean flux density to match the mean of the "starless" 
group. We flnd 52%, 64%, and 91% for groups 1 — 3, re- 
spectively. This is an order of magnitude higher than the 
mean free-free flux density contribution we find with the 
RMS 6 cm observations. Therefore, it is highly unlikely 
that the trend of increasing fiux density is due to free- 
free contamination. However, free-free emission could 
contribute significantly to the flux densities of individual 
sources. 

The observed increase in flux density with group could 
also be a reflection of an increasing dust temperature. 
While an increase in does explain some of the ob- 
served trend in flux density, it cannot explain the or- 
der of magnitude increase in mean flux density between 
groups and 3. We discuss the effects of To in Section 
below. 

A SIMBAD search of the 22 "starless" sources in the 
inner Galaxy with S^> 10 Jy revealed that only 4 of 
the brightest BGPS sources are possibly truly starless, 
while the remaining 18 brightest BGPS sources in the 
inner Galaxy had a combination of radio emission and 
masers as indicators of active star formation. This ac- 
centuates that while these BGPS sources are "starless" 
in the sense that they do not contain any mid-IR sources 
from GLIMPSE or RMS catalogs, they are not necessar- 
ily truly starless sources. Further comparison to other 
surveys will be needed in order to find a truly starless 
sample of BGPS sources. 

Figure [6] shows the number of BGPS sources as a func- 
tion of the concentration parameter in each star forma- 
tion probability group. The concentration parameter is 
defined as the ratio of the elliptical radius containing 90% 
of the total source flux to the elliptical radius contain- 
ing 50% of the total source flux, where the elliptical ra- 
dius is given hy R — {RmajRminY^^ ■ By this definition, 
compact sources will have a low concentration parameter 
while sources with large skirts of emission will have a high 
concentration parameter. The "starless" sources have a 
lower mean concentration parameter, and the mean value 
increases with probability of star formation activity. The 
majority of the most extended BGPS sources fall within 
probability Group 3. These statistics are likely affected 
by chance alignments since more extended sources are 
more likely to be identified by chance alignments with the 
mid-IR sources. Since we cannot identify which BGPS- 
mid-IR source associations are due to chance alignments. 
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Fig. 6. — Number of sources as a function of BGPS concentration parameter for the different probability groups of BGPS sources. Only 
the 5067 BGPS sources in the region where all catalogs overlap are included. The concentration parameter is given by the ratio of the 
elliptical radius containing 90% of the millimeter emission to the elliptical radius including 50% of the emission. A smaller concentration 
parameter means the source is more compact, while sources with large skirts of emission will have a higher concentration parameter. 
Note the vertical scale is different for each probability group. The range of confidence values for each probability group are given in the 
upper-right corner of each panel. 



we cannot assess the extent of the effect of chance ahgn- 
ments on this trend. However, we do note that the group 
3 sources are unhkely to include many chance ahgnments 
since we expect a very low fraction of matched EGO and 
RMS sources to be the result of chance alignments. 

Another possible bias which could affect the statistics 
of each probability group is distance. Due to the limited 
sensitivity of the mid-IR catalogs, it is possible that more 
distant BGPS sources harbor low-mass YSOs which are 
not present in the mid-IR catalogs. The starless group 
could include a large number of distant BGPS sources 
which are actually forming stars. 

4.2. Physical Properties 

We create a distance subsample of BGPS sources 
within the region where all surveys overlap that have 
distance estimates from a molecular line study of BGPS 
sources (Schlingman et al. 2011) and VLBI parallax mea- 
surements (7 BGPS sources; Bartkiewicz et al. 2008; 
Brunthaler et al. 2009; Xu et al. 2009). Schlingman et 
al. used the VLBI parallax based model of Galactic ro- 
tation of Reid et al. (2009) and were able to assign the 
near kinematic distance to sources based on the asso- 
ciation with a GLIMPSE IRDC (Peretto et al. 2009). 
Schlingman et al. additionally adopted distances from 
the well-studied sample of massive star-forming regions 
of Shirley et al. (2003). The "distance subsample" pre- 
sented in this section consists of 280 BGPS sources which 
are located in the region where all the catalogs overlap. 
Since the dominant number of distance determinations 
are based on the association of a BGPS source with an 
IRDC, the distance subsample is significantly biased to 
the near kinematic distance in the inner Galaxy. Full re- 



sults of the molecular line study are formally presented 
in Schlingman et al. (2011). 

We calculate properties such as mass, volume-averaged 
density, column density, etc., for the 280 BGPS sources 
with kinematic distances. The distance subsample sam- 
ples the different probability groups well with 55, 67, 97, 
and 61 from groups 0-3, respectively. Figure [7] plots a 
histogram of the distances as a function of probability of 
star formation activity. The distances range from 0.21 to 
7.0 kpc with no obvious trends with star formation prob- 
ability (see Table 131) ■ Since the BGPS will be sensitive to 
different structures over such a large range in distances, 
we define a subsample of sources with 2 kpc< d < 3 
kpc (hereafter refered to as the 2 kpc subsample) in or- 
der to separate trends from the ambiguity introduced by 
a large range of distances. We have chosen this range 
of distances because Dunham et al. (2010) find that at 
distances greater than approximately 2 kpc the BGPS 
is sensitive to structures most closely related to clumps. 
Thus, we expect the 2 kpc subsample to have properties 
in the typical range of clumps (e.g. Table 1 of Bergin 
& TafaUa 2007). We find 6, 10, 18, and 8 sources with 
distances of 2 — 3 kpc for probability groups 0-3, respec- 
tively. 

Figure [8] plots the number of objects versus physi- 
cal radius for each probability group. The gray his- 
tograms show the distance subsample, while the striped 
histograms mark the 2 kpc subsample. The elliptical 
radius was determined by the BGPS source extraction 
software based on the emission weighted moments of the 
BGPS maps (Rosolowsky et al. 2010). Radii of unre- 
solved sources have been set to an upper limit of half the 
beam size, 16.5". The physical radii of source in the dis- 
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Fig. 7. — Number of sources as a function of distance for each of tiic different probability groups of BGPS sources which have distance 
measurements and are located in the region where all catalogs overlap. The striped histogram represents the 2 kpc subsample. Note the 
vertical scale is different for each probability group. The range of confidence values for each probability group are given in the upper-right 
corner of each panel. 
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Fig. 8. — Number of sources as a function of physical radius for each of the probability groups of BGPS sources which have distance 
measurements and are located in the region where all catalogs overlap. The striped histogram represents the 2 kpc subsample. The range 
of confidence values for each probability group are given in the upper-right corner of each panel. 
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tance subsample range from 0.07 to 3.8 pc spanning the 
range of sizes expected for the entire hierarchical struc- 
ture in GMCs including clouds, clumps, and cores (radius 
= 1 — 7.5 pc, 0.15 — 1.5 pc and 0.01 — 0.1 pc, respectively; 
Bergin & Tafalla 2007). The radii of the 2 kpc subsample 
range from 0.18 — 1.4 pc, and span the range of typical 
clump radii (see Table [3]) . The physical radius increases 
with star formation probability group from 0.95 pc for 
the "starless" sources to 1.3 pc for group 3 (see Table [3]). 

We calculate the isothermal dust mass, M^so, assuming 
a single dust temperature, Tu, as 

M- 

-""ISO ,rTi \ 

= 13.1Mo (^,13.0 K/20.0A'_ij' 

(1) 

where S^, is the flux density, D is the distance, Ki, is the 
dust opacity per gram of gas and dust including a gas- 
to-dust ratio of 100 {k^, = 0.0114 cm^ g-i at 271.1 GHz; 
logarithmically interpolated from 0H5 dust from column 
5 of Table 1 in Ossenkopf & Henning 1994), and is 
the Planck function evaluated at Tjj. The right hand 
side of the equation is evaluated for T£)= 20 K based on 
NH3 observations of BGPS sources in the Gemini OBI 
molecular cloud where a mean gas kinetic temperature 
of 20 K was found (Dunham et al. 2010). 

Although 20 K is a reasonable value of Tu to assume 
if no other information is available, we have both gas ki- 
netic temperature and dust temperature measurements 
from other observations. M. K. Dunham et al. (in prepa- 
ration) have surveyed 716 BGPS sources in the inner 
Galaxy in NH3(1,1), (2,2), and (3,3) and have provided 
a measurement of Tk towards the 408 BGPS sources de- 
tected in NH3. They find (Tk) of 14.0±3.2 K, 14.8±3.9 
K, 16.1±4.3 K, and 22.0±5.9 K based on 130, 143, 102, 
and 33 BGPS sources from groups — 3, respectively. For 
this work we assume the mean gas kinetic temperatures 
from the NH3 study for each group. By assuming the gas 
kinetic temperture is equivalent to the dust temperature, 
we are assuming that the density is high enough to bring 
the temperatures into equilibrium via collisions. This 
assumption holds in the dense regions where the dust 
effectively shields the gas from the interstellar radiation 
field (ISRF) and the gas is heated via collisions with the 
dust (e.g. Goldreich & Kwan 1974; Evans et al. 2001). 

The number of BGPS sources of each Miso are plotted 
for each probability group in Figure[Hl where the gray his- 
togram represents the entire distance subsample and the 
striped histogram denotes the 2 kpc subsample. There 
is a clear trend of increasing mass with increasing prob- 
ability of star formation activity; the median masses in 
groups — 3 are 340, 690, 680, and 1160 M©, respectively. 
Similarly, we find a median mass of 52, 94, 170, 710 Mq 
for groups — 3 of the 2 kpc subsample, respectively. The 
median mass is a more accurate representation of each 
group than the mean mass due to the few extremely high 
mass (10** Mq) BGPS sources in groups 0, 2, and 3 (see 
Table O. 

There are two effects that could artificially produce 
the observed trend in mass: contamination by free-free 
emission and the assumed T^. As discussed previously 



in section 14.11 the free-free emission is not a significant 
contribution to the flux densities of most BGPS sources 
and is unlikely to cause the trend of increasing mass. 

Here we explore what effects different dust tempera- 
tures would have on the observed trend in mass. We 
calculate the temperature required for the median mass 
of each group to equal the median mass of group 1 (where 
we assume Tu= 14.8 K). We have chosen to use the me- 
dian mass of group 1 as the standard since it is very 
similar to the median for group 2 (690 and 680 M©, re- 
spectively). Dust temperatures of 9.2 K, 15.9 K, and 
33.2 K would bring the median of groups 0, 2, and 3 
into alignment with group 1, respectively. While these 
temperatures are reasonable, the values for groups and 
3 lie outside of the standard deviation of the mean gas 
kinetic temperatures derived from the NH3 observations 
discussed above. The NH3 observations are sensitive to 
gas kinetic temperatures well above 33 K, so the lower 
mean temperature observed toward group 3 sources is 
not the result of an observational bias. 

Thus, the trend to higher masses with probability of 
star formation activity is likely real. Miso ranges from 
2.4 Mq to 3.6 X 10^ Mq for the distance subsample, and 
19 Mq to 2.5 X 10^ M© for the 2 kpc subsample. Again, 
the distance subsample spans the range of masses for 
cores, clumps, and clouds (0.5 — 5 M©, 50 — 500 M©, 
and 10^ - 10" M©, respectively; Bergin & Tafalla 2007) 
while the 2 kpc subsample encompasses the mass range 
for clumps without extending into the range for clouds 
or cores. 

From Miso we calculate the mean particle density, 
n, for each source in the distance subsample via n = 
SMiso/iTrR^ , where R is the source radius determined 
by the BGPS source extraction algorithm. This density 
calculation assumes a constant density and the effect is 
a volume-averaged density that is more reflective of the 
lower density outer regions of the BGPS source than the 
higher density regions. Figure [TOl shows histograms of n 
for each probability group. All groups span almost the 
same range from 320 to 5.6 x 10"* cm"'^. However, the 
BGPS sources in groups and 3 are, on average, higher 
density than those in groups 1 and 2 (see Table |3]) . The 
"starless" sources are less massive but also smaller than 
BGPS sources in groups 1 — 3 and therefore have densities 
similar to those of group 3. The 2 kpc subsample only 
extends to about 10 cm^"^ on the low density end of each 
distribution, but does extend to 10"* cm^'^ on the high 
density end, even in the "starless" category. Clouds typ- 
ically range from 50 — 500 cm~^, clumps from 10^ — lO" 
cm""^, and cores from 10"*^ — 10^ cm~'^ (Bergin & Tafalla 
2007). Overall, the distance subsample spans the entire 
range of densities, but the 2 kpc subsample only spans 
the density ranges for clumps and cores. 

The column density of H2 is given by N'^f — 

Aliso/'i'T^J.mHRobj^ where = 2.37 and rriH is the mass 
of hydrogen. Figure [TT] shows the number of sources as 
a function of N^^^^ for each probability group in the dis- 
tance subsample (gray histogram) as well as for the 2 kpc 
subsample (striped histogram). There is an increase in 
mean N"^^^ between groups 0—2 and Group 3. Groups 

— 2 have mean N?^^ values near 1.2 x 10^^ cm~^ while 
Group 3 has a mean H2 column density of 2.4 x 10^^ 
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Fig. 9. — Number of sources as a function of isothermal mass for each of the probability groups of BGPS sources which have distance 
measurements and are located in the region where all catalogs overlap. The striped histogram represents the 2 kpc subsample. Note the 
vertical scale is different for each probability group. The range of confidence values for each probability group are given in the upper-right 
corner of each panel. 




log(n) [cm ^] 



Fig. 10. — Number of sources as a function of mean volume-averaged density for the subset of each BGPS source probability group 
which have distance measurements and are located in the region where all catalogs overlap. The striped histogram represents the 2 kpc 
subsample. Note the vertical scale is different for each probability group. The range of confidence values for each probability group are 
given in the upper-right corner of each panel. 
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Fig. 11. — Number of sources as a function of H2 column density for the subset of each BGPS source probability group which have 
distance measurements and are located in the region where all surveys overlap. The striped histogram represents the 2 kpc subsample. 
The dotted line denotes 1 g cm^^, which is the theoretical minimum surface density needed to prevent fragmentation and allow a massive 
star to form (Krumholz & McKee 2008). The dashed line denotes 122.5 Mq pc~^ which is the threshold column density above which 
star formation occurs that has been observed in nearby low-mass star-forming regions (Lada ct al. 2010; Heiderman et al. 2010). Note the 
vertical scale is different for each probability group. The range of confidence values for each probability group are given in the upper-right 
corner of each panel. 




Fig. 12. — Number of sources as a function of H2 column density per beam for the subset of each BGPS source probability group for all 
BGPS sources within the region where all catalogs overlap (10° < £ < 65.5°). The dotted line denotes 1 g cm~^, which is the theoretical 
minimum surface density needed to prevent fragmentation and allow a massive star to form (Krumholz & McKee 2008). The dashed 
line denotes 122.5 Mq pc~^ which is the threshold column density above which star formation occurs that has been observed in nearby 
low-mass star-forming regions (Lada et al. 2010; Heiderman et al. 2010). Note the vertical scale is different for each probability group. The 
range of confidence values for each probability group are given in the upper-right corner of each panel. 
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cni^^(see Table [3]). The dashed hne denotes the thresh- 
old column density of about 123 Mq pc~^ which has been 
observed toward nearby, low-mass star-forming regions 
(Lada et al. 2010; Heiderman et al. 2010). The dotted 
line denotes the theoretical minimum surface density, 1 
g cm~^, required to prevent fragmentation and allow a 
massive star to form (Krumholz & McKee 2008). Similar 
to the volume density, N'^^ is also an average quantity. 
There are likely higher surface density regions within the 
BGPS source where the YSOs are forming, but the BGPS 
is sensitive to the lower density outer regions of the clump 
(Dunham et al. 2010). This is consistent with the pic- 
ture that RMS and EGO sources are high-mass YSOs 
even though the average surface density is not above the 
1 g cm^^ requirement to prevent fragmentation. 

Alternatively, we can calculate the H2 column density 
per beam for all BGPS sources within the area where 
all catalogs overlap since this quantity is independent of 
distance. The H2 column density per beam is given by 

jybeam _ (^0 ) ^2 <j 

where Si,{AQ") is the 1.1 mm flux density within an aper- 
ature with a diameter of 40", ^beam is the solid angle of 
the beam, = 2.8, k.^, is the dust opacity and B^(Tu) 
is the Planck function evaluated at the dust temperature, 
Td, as described above. We use 5*1,(40") as a measure of 
the flux within a beam since a top-hat function with a 
20" radius has the same solid angle as a Gaussian beam 
with a FWHM of 33". We have applied an aperture cor- 
rection of 1.46±0.05 to 5*1/(40") to account for power that 
falls outside of the 40" aperture due to the sidelobes of 
the CSO beam (Aguirre et al. 2011). We calculate this 
property for all BGPS sources within the region where all 
surveys overlap since N^^"^ is independent of distance. 
Figure [12] shows the distribution of H2 column density 
per beam for the BGPS sources in the overlap region 
for each probability group. As above, the dotted line 
denotes the column density required to prevent fragmen- 
tation (Krumholz & McKee 2008) and the dashed line 
denotes the threshold for effecient star formation (Lada 
et al. 2010; Heiderman et al. 2010). The mean column 
density per beam increases with probability group with 
^jY&eam^^ 0.95x10^2, 1.4x10^2, l.GxlO^^, and 4.7x10^2 

cm~2 for groups — 3, respectively. The median jy'^^"'- 
of group (7.1 X 10^^ cm~2) is similar to the threshold for 
efficient star formation (6.5 x 10^^ cm~^; dashed line in 
Figure [T^ Lada et al. 2010; Heiderman et al. 2010) such 
that roughly half of the "starless" BGPS sources have 
a peak column density high enough to efficiently form 
stars. It may be that these sources appear "starless" in 
this study because of the sensitivity limits of the mid-IR 
catalogs. For example, low-mass YSOs within a BGPS 
source at 5 kpc would not be detected by the GLIMPSE 
surveys. Further comparisons with other tracers of star 
formation are necessary to determine which sources are 
truly starless. 

OveraU, we find that 5„(int), R, Mi,o, N^f^"", and 
N"^^^ are all higher in BGPS sources actively forming 
stars. Sources identified as actively forming stars are 
typically extended with large skirts of emission, while the 
"starless" sources are more compact, although this trend 



is likely affected by chance alignments since more ex- 
tended BGPS sources are more likely to be matched with 
a mid-IR source by chance. BGPS sources with EGO 
or RMS sources (group 3) are more massive and have 
higher volume-averaged densities than BGPS sources 
with GLIMPSE sources. As expected due to the large 
range of distances, these physical properties show that 
the BGPS sources span the entire hierarchical structure 
in GMCs (clouds, clumps, and cores). Based on the phys- 
ical properties presented, we confirm that the 2 kpc sub- 
sample is mostly comprised of clumps. 

5. SUMMARY 

We have cross-matched the Bolocam Galactic Plane 
Survey (BGPS) catalog with available mid-infrared 
Galactic plane catalogs to determine which BGPS 
sources show signs of active star formation. The mid- 
infrared catalogs include the GLIMPSE Red Source cat- 
alog (R08; Robitaille et al. 2008), the Extended Green 
Object catalog (EGOs; Cyganowski et al. 2008), and the 
Red MSX Source catalog (RMS; Hoare et al. 2004). We 
have also created an additional GLIMPSE Red Source 
catalog with a less stringent color selection than was used 
by R08, namely [4.5] - [8.0] > 0.75. 

OveraU, 3,712 BGPS sources (44%) include at least 
one mid-IR source. Within the area where all surveys 
overlap (10° < £ < 65°) 2,457 of 5,067 (48%) BPGS 
sources include at least one mid-IR source. If we ac- 
count for the number of expected chance alignments be- 
tween the mid-IR and BGPS catalogs, we estimate the 
total number of BGPS sources showing mid-IR signs of 
active star formation to be 1,472 over the whole sur- 
vey (17.6%) and 1,035 within the area where all surveys 
overlap (20.4%). These low percentages are conserva- 
tive lower-limits, while the true fraction of BGPS sources 
forming stars is likely higher. 

We additionally divide the BGPS sources into four 
groups representing the probability of star formation ac- 
tivity based on the highest confidence parameter mid- 
IR source matched to each BGPS source in order to 
study physical properties as a function of star forma- 
tion probability. The sources with the highest probabil- 
ity of star formation activity are BGPS sources with an 
EGO or young RMS source (90% confidence). The low- 
est probability group corresponds to sources which were 
not matched with any mid-IR sources and are considered 
to be "starless" in this picture. 

The mean BGPS flux increases with probability of star 
formation activity, with the BGPS sources matched to 
EGOs and young RMS sources statistically having the 
highest fluxes. The distribution of "starless" sources 
has the lowest mean flux but does include a tail up to 
^30 Jy. While BGPS sources without mid-IR associa- 
tions are designated "starless" here, a literature search 
shows that these bright sources do actually show signs 
of active star formation in other tracers, such as masers. 
X-ray emission and radio emission, which demonstrates 
the importance of including other star formation tracers 
when defining a truly starless sample. For example, these 
"starless" sources may harbor young, highly embedded 
YSOs which only exhibit emission at wavelengths long- 
ward of 24 /Ltm. A future paper will compare the BGPS 
catalog with the Spitzer MIPSGAL catalog of 24, 70, and 
160 /Ltm sources. Additionally, this literature search sug- 
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gests that the fraction of BGPS sources with active star 
formation is hkely higher than our conservative lower- 
hmit of 18%. 

The most compact BGPS sources are typically found 
to be "starless", while BGPS sources with the highest 
probability of active star formation are typically more 
extended with large skirts of emission. 

We define a sample of 280 sources with kinematic dis- 
tance determinations ranging from 0.21 to 7.0 kpc, which 
are strongly biased to the near distance in the inner 
Galaxy. Since BGPS emission arising from 7 kpc will 
trace a different average density of gas than that arising 
from 1 kpc, we define a subsample of sources with dis- 
tances between 2 and 3 kpc. We see a trend of increasing 
mass with increasing probability of star formation activ- 
ity, both in the distance sample and the 2 kpc subsample. 
Similarly, there is a trend in H2 column density, with the 
highest mean column density found in the group with the 
highest probability of star formation. 

The physical properties presented show that the BGPS 
sources span the entire hierarchical structure in GMCs 
(clouds, clumps, and cores), as expected due to the large 
range of distances. We also confirm that the 2 kpc 
subsample is mostly comprised of clumps. On aver- 
age, BGPS sources matched with EGO or young RMS 
sources are more massive and denser than BGPS sources 
matched with GLIMPSE red sources or not matched with 
any mid-IR sources. 
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TABLE 1 
Matching Statistics 
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33.5 


8.2±0.4 


7360 
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EGOs 


84 
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0.13±0.05 
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RMS total 
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8.4±1.2 


6062 
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0.78±0.11 


88.6 


Unclassified / Other 


26 


13 


50.0 


7.2±4.7 


6062 


13 
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0.029±0.019 


86.5 


Evolved Star 
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Young/old star? 
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HII Region 
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Total 
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44.4 
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TABLE 2 

MiD-IR CONTENT OP THE BGPS SOURCES 
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-28.99171 


0.01896 


-0.13780 


82.1 











3 


GOOO. 010+00.157 


266.26410 


-28.83834 


0.01913 


0.15628 


81.4 








7 


3 


GOOO.016-00.017 


266.42551 


-28.93289 


0.01216 


-0.01361 


93.2 











6 


GOOO.018-00.431 


266.83618 


-29.14625 


0.01632 


-0.43125 


< 16.5 











2 


GOOO.020+00.03.3 


266.38130 


-28.90146 


0.01881 


0.03579 


77.2 











7 


GOOO. 020-00. 051 


266.468.30 


-28.94505 


0.02127 


-0.05191 


62.3 








1 


1 


GOOO.022+00.251 


266.17683 


-28.78469 


0.02494 


0.24951 


< 16.5 








1 


1 


GOOO.034-00.437 


266.85490 


-29.13193 


0.03704 


-0.43782 


< 16.5 















Note. — The full table is available in the online Journal. 
^ The BGPS source name is based on the Galactic coordinates of the peak of the 1.1 mm emission. 

The coordinates of the geometric centroid of the 1.1 mm emission. 
^ For sources unresolved compared to the beam, the radius is set to an upper-limit equal to half the beam size, 
16.5". 

^ Only young RMS types ("HII Region", "HII/YSO", "YSO", and "Young/old star?") are included here. 
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TABLE 3 

Physical Property Statistics of BGPS sources in Region Where All 
Catalogs Overlap 



Property 


1 1 {JUaUlLi ty 






ci" Q Ti a v/H 
btctllUdil U 






(units) 


Group 


miiuntium 


mean 


deviation 


median 


maximum 







0.06 


0.71 


1.4 


0.41 


33.1 


(Jy) 


1 


0.06 


1.5 


1.9 


0.89 


27.7 


2 


0.08 


2.0 


2.5 


1.2 


25.3 




3 


0.13 


8.0 


18.0 


3.4 


146.8 


(Jonccntration 





1.3 


1.6 


0.11 


1.6 


2.1 


Parameter 


1 


1.3 


1.7 


0.13 


1.6 


2.6 




2 


1.4 


1.7 


0.13 


1.7 


2.3 




3 


1.5 


1.8 


0.20 


1.8 


2.7 


j\jbeam 





0.21 


0.95 


1.2 


0.71 


23.3 


\ 


0.20 


1.4 


1.7 


0.95 


27.0 


2 


0.22 


1.6 


1.6 


1.1 


19.3 




3 


0.20 


4.7 


10.4 


2.1 


99.3 






Distance Subsampl 


e 






Distance 





1.04 


4.3 


1.2 


4.4 


6.7 


(kpc) 


1 


0.21 


4.2 


1.2 


4.5 


6.4 




2 


1.9 


4.1 


1.2 


3.8 


7.0 




3 


1.3 


4.4 


1.4 


1.6 


7.0 


Kadius 





0.12 


0.95 


0.55 


0.86 


2.0 


(pc) 


1 


0.07 


1.1 


0.63 


1.0 


3.4 


2 


0.18 


1.2 


0.75 


1.1 


3.5 




3 


0.35 


1.3 


0.72 


1.1 


3.8 







19 


1050 


1720 


340 


1.0x10* 


(Mo) 


1 


2.4 


1040 


1130 


690 


5380 


2 


30 


1450 


2110 


680 


1.2x10" 




3 


90 


3230 


5450 


1160 


3.6x10-* 


n 





0.73 


7.1 


11.0 


3.3 


56.1 


(103 cm-3) 


1 


0.32 


6.7 


9.4 


2.5 


50.2 


2 


0.41 


5.8 


10.7 


2.5 


85.7 




3 


0.60 


7.2 


9.7 


3.2 


45.8 


Ho 

(lO^i cm-2) 





3.1 


12.8 


9.7 


9.6 


56.3 


1 


3.7 


12.8 


9.9 


9.7 


47.8 


2 


3.2 


12.4 


10.8 


9.2 


82.2 




3 


3.4 


24.2 


31.3 


12.6 


154.0 


2 kpc Subsample 


Kadius 





0.19 


0.49 


0.29 


0.51 


0.90 


(pc) 


1 


0.19 


0.65 


0.41 


0.49 


1.4 


2 


0.18 


0.67 


0.27 


0.65 


1.1 




3 


0.41 


0.74 


0.23 


0.73 


1.1 


JWiso('Jb) 





37 


79 


56 


52 


179 


(Mo) 


1 


26 


390 


767 


94 


2520 


2 


19 


192 


130 


170 


546 




3 


90 


792 


609 


710 


2120 


n 





1.0 


8.2 


11.3 


3.6 


29.5 


(103 cm-3) 


1 


0.85 


5.3 


5.8 


3.7 


16.3 


2 


0.91 


4.0 


4.2 


2.3 


14.2 




3 


2.0 


9.1 


6.9 


8.1 


22.3 


(10^1 cm-2) 





2.7 


7.2 


6.8 


4.7 


19.9 


1 


2.4 


7.9 


6.9 


6.1 


25.4 


2 


3.1 


6.4 


3.5 


5.0 


14.5 




3 


6.3 


21.0 


13.7 


20.4 


44.6 



